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We study the equation of state in two-flavor QCD at finite temperature and density. Simula- 
tions are made with the RG-improved gluon action and the clover-improved Wilson quark action. 
Along the lines of constant physics for mps/my = 0.65 and 0.80, we compute the derivatives of 
the quark determinant with respect to the quark chemical potential fj, q up to the fourth order at 
flq = 0. We adopt several improvement techniques in the evaluation. We study thermodynamic 
quantities and quark number susceptibilities at finite ji q using these derivatives. We find enhance- 
ment of the quark number susceptibility at finite \i q , in accordance with previous observations 
using staggered-type quarks. This suggests the existence of a nearby critical point. 
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Figure 1: Phase diagram and lines of constant physics at /J. = for two-flavor QCD with improved Wilson 
quarks and Iwasaki glue at N t = 4 [g, |[] . 

1. Introduction 

Finite density QCD has been studied on the lattice mainly using staggered-type quarks [jl]]. 
However, because the expected 0(4) universality of the deconfming transition in two-flavor QCD 
has not been confirmed with staggered-type quarks, the results may contain sizable lattice artifacts. 
Therefore, we need to crosscheck the results with different lattice actions. We study it with Wilson- 
type quarks, with which the 0(4) scaling has been confirmed [Q, ||]. 

Because Wilson-type quarks are numerically more intensive, we have to adopt/develop several 
improvement techniques. We apply a hybrid method of Taylor expansion and spectral reweighting, 
and develop a couple of improvement tricks. 

2. Formulation 

We extend the study of two-flavor QCD by the CP-PACS Collaboration at vanishing chemical 
potential \i q = [|3|] to finite densities. Preliminary reports of this study have been presented at 
Lattice 2006 and 2007 conferences @]. 

We adopt the clover-improved Wilson quark action and the RG-improved Iwasaki gauge action 
defined by S = S g + S q with 

S g = -PZ\ C ° I ^(xj + ci t ^v X2 M], S q =^4M x ,qf, (2.1) 
x \ p.<v;ji,v=\ n=£v;ii,v=l j /=1 x,y 

where j3 = 6/g 2 , c\ = —0.331, Co = 1 — 8ci, and 

M x , y = d^-K^iil-YdUxjS^ + il + YdUlA^} 

i=\ 

-K{e»{\ - Y,)U xA 8 x+ly + e^{\ + 74)^4^+4} " I <W- ( 2 - 2 ) 

where /J. = \i q a and F^ v is the lattice field strength in terms of the standard clover-shaped com- 
bination of gauge links. For the clover coefficient csw> we adopt a mean field value using W lxl 
calculated in the one-loop perturbation theory: c sw = (W lxl )~ 3/4 = (1 — 0.8412j8 ^) 3/4 [gj. 
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The lattice size is xN t = N s - lte = 16 3 x 4. We carry out simulations along the lines of con- 
stant physics (LOCs) for mps/my = 0.65 and 0.80 at jj, q = ^ shown in FigJI[ For definiteness, 
we refer the pseudocritical temperature T pc at /I = as Tq in the followings. Runs are carried out 
in the range j3 = 1.50-2.40 at thirteen values of T /Tq ~ 0.82-4.0 for mp S /m v = 0.65 and twelve 
values of T /Tq ~ 0.76-3.0 for nips/my = 0.80. The number of trajectories for each run after 
thermalization is 5000-6000. We measure physical quantities at every 10 trajectories. 



3. Taylor expansion up to 0(/i 4 ) 

We expand the grand canonical potential ft) = p/T 4 in terms of \l q /T at jj, q = 0: 



CO 



1 N? d n ln3f 



n\Nf d{n q /Ty 



(3.1) 



where V = (N s a) 3 is the system volume and 3? is the partition function. We expand up to (pL q /T) 
The relevant coefficients are given by 



N t 

2NT 



1 



? 2 , C 4 



4W s 3 iV ( 



2 2 ), 



(3.2) 



^2 = (@2) + <^i 2 > , < = {%) +4(^ 3 ^i) + 3 + 6 (^ 2 ^ 2 > + {9f) , 



where !3 n = Nf ij^j IndetM is calculated as 



— , S> 2 =Nf 



tr I M 



d 2 M 



tr M" 



dM ,dM 

M^ 1 

d}X d}X 



etc. 



We calculate similar coefficients for the expansion in terms of the isospin chemical potential jij 
too. 

We evaluate the traces with the random noise method. In the evaluation, we apply the fol- 
lowing two improvements: (i) Because the elements off-diagonal in color and spin indices are not 
suppressed by \x — y\ with Wilson-type quarks, the number of the same-magnitude off-diagonal el- 
ements in the quark matrix is 1 1 times larger than the diagonal one. This is different from the case 
of staggered-type quarks, in which off-diagonal elements in spin indices are slightly suppressed by 
the spatial offset. Because a large number of noises is required to pick up a signal from data with 
S/N = 1/1 1, we decide not to apply the noise method for color and spin indices and generate noise 
vectors only for spatial indices, i.e. we repeat the inversion of M for each color and spin indices, 
(ii) We find that the dominant errors are from terms containing Therefore, we adopt 10-40 
times more noise vectors for &\. while we generate only 10 noise vectors for other traces. 

Results for the pressure difference Ap = p(pL q ) — p(0), quark number density (n q ), and quark 
number as well as isospin susceptibilities (% q and Xi, respectively) are shown in Fig^for raps/ray = 
0.65. We observe larger enhancement in the quark number fluctuations near Tq with increasing \i q , 
however such an enhancement around Tq is not shown in the isospin fluctuations. Because the 
statistical errors in n/T 3 and % q jT 2 are large, further studies are needed increasing the statistics for 
more precise arguments with this approach. 



3 



EOS at finite density in two-flavor QCD 



K. Kanaya 




Figure 2: Results of the Taylor expansion method up to the fourth order of chemical potentials for mps / my = 
0.65. Tq is the pseudocritical temperature at jj. — for the same LOC. 



4. Improvement: A Hybrid Method 



In the previous section, the Taylor expansion is limited up to the fourth order of the chemical 
potential, because the evaluation of & n with n > 4 is compuationally demanding. To improve the 
calculation, we need to estimate c„ at larger n. Here, we note that @„ = at n > 4 for the free quark 
case. Therefore, at high temperatures, we may approximate Ss n = for n > 4 in the evaluation of 
c n at n > 4. This approximation corresponds to a hybrid reweighting method in which the grand 
canonical potential is approximated by a truncated Taylor expansion 



with N„ 



with 



<o(T,n) 



4. Here, 



1 



VT 3 



injr(r,o) + 



^In^exp 



n=l 



(M=0) 



>(u=0) 1S the average over configurations at /I = 0. We then have 

3T(T,n) w iT(r,0) ( gf (M)£,!e(Al) ) (/J _ 0) 



1 



AT f Im[lndetM(ju)] w Im^i/x + ^ylm^ 3i u 3 , 



M'Re 



ln / detM( M ) \ 
VdetM(O) / 



^Re^ 2A t 2 + ^Re^4At 4 - 



(4.1) 

(4.2) 

(4.3) 
(4.4) 



This kind of hybrid method was first tested in [^] with two flavors of staggered quarks. In that 
study, it turned out to be difficult to control statistical errors at \i q larger than 0(T). This is due to 
the sign problem at finite \i q . 
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Figure 3: Typical 9 -distribution of two-flavor QCD at ji q = for wps/my = 0.65. 
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Figure 4: Binder cumulant test of 9 -distribution for raps /my = 0.65 (left) and 0.80 (right). Data labeled as 
"no sub." are the results of approximate estimation of 9 2 and 9 4 without subtracting the contribution of the 
same noise vectors. Data without "no sub." are the correct estimates. 



In a previous study with staggered quarks, SE noted that the 6 -distribution at \x = is well 
described by a Gaussian form, and showed that this fact can be used to carry out the 6 -averaging 
with small errors [0|. Furthermore, we can argue that, because 6 = &(n) at small jj,, deviations 
from the Gaussian distribution do not affect the terms up to }X A We find that our data are 

also well Gaussian (see Fig.|| for typical examples). A convenient test of Gaussian distributions 
is provided by the Binder cumulant B® = (0 4 ) / (0 2 ) 2 - If the distribution is Gaussian, we expect 
B4 = 3. In Fig-H we show that all of our distributions are approximately consistent with B% = 3. 

With Gaussian 0-distribution, we can carry out the ^-averaging, and the task is reduced to 
evaluate 

3r(T,n)*3r(T,0)( < ?e-tt e2 )A i , (4.5) 
\ / (ju=0) 

with (0 2 ) F the mean d 2 for given fixed F . 

The remaining F-averaging is again challenging because the factor e F ^ can easily shift the 
central contribution for the average to a statistically poor region of F. At small /i, this problem can 
be largely resolved by shifting j8 with /I such that the fluctuation in e F ^ is compensated by that in 
the gauge action, because F is sensitively correlated with the gauge action S g . More concretely, we 
adopt the reweighting method to shift j8 and calculate an optimal j3 by minimizing the fluctuation 

/ (e F e-^ e2 )Fe 6N ^-^ p - (^e~^)' g^tf-W^X 
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Figure 5: Results of EOS using the hybrid method for mps/my = 0.65. 
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Figure 6: Results of EOS using the hybrid method for raps /my = 0.80. 



where P = S g / (6A^ s i te j8) is the generalized plaquette. The shifts in j8 turn out to be less than about 
0.03 in our study. Since these are negligible in Fig.[l|, we disregard the resulting small deviation 
from the line of constant physics, and simply translate the shifts in j8 to shifts in T for the final 
plots. 

Our results of EOS with these improvements are summarized in Figs.|| and |6| for mps / my = 
0.65 and 0.80, respectively. Here, we calculate the quark number density and its susceptibility 
by numerical differentiations of the grand canonical potential using the following thermodynamic 
formulae: 

n q _Nf apnjT) Xq _N? d 2 (\n$f) 

We find that, in spite of the fact that the simulations at different T /To are independent, the resulting 
EOS shown in these figures has smooth T- and ji q -dependence which is in accordance with theoret- 
ical expectations. Therefore, we think that the assumption & n = at n > 4 is well satisfied down to 
T close to the transition temperature. Furthermore, we find that the statistical fluctuations are much 
reduced over the results of the previous section (Fig|2]). This is due to the improvements including 
the Gaussian method for ^-averaging and the j3 -reweighting for F-averaging. The hybrid method 
with improvements provides us with a tractable way to calculate EOS with Wilson-type quarks. 

Similar to the case of the previous section, the quark number susceptibility for raps/ray = 0.65 
show enhancement near To with increasing pL q , although the statistical errors are not quite small yet. 
On the other hand, the quark number susceptibility for mps/my = 0.80 shown in Figj^do not show 
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rapid enhancement with \x q . This may be in part explained by the expectation that the critical 
point locates at larger \i q because the quark mass is larger than that for mps/ray = 0.65. The mild 
enhancement shown in Figj6| may be suggesting that, at this quark mass, the critical point does 
not locate inside the applicability range in \l q with our &{H q ) calculation. Further studies with 
increased statistics around 7b are needed for more definite conclusions. 



5. Conclusions 



We have carried out the first calculation of the equation of state at non-zero densities with two 
flavors of improved Wilson quarks. Statistical fluctuations of physical observables at finite density 
are much severer with Wilson-type quarks than with staggered-type quarks. To tame the problem, 
we combined and developed several improvement techniques. 

With these improvements, we found that the peak height of the quark number fluctuation at the 
pseudo-critical temperature increases as \i q increases. In contrast, isospin susceptibilities show no 
sharp peaks at the pseudo-critical temperature. These results agree with previous observations by 
the Bielefeld-Swansea Collaboration using staggered-type quarks, and suggest that a critical point 
exists at finite pL q , which is expected to locate at the end point of a first order transition line between 
confining and deconfining phases in the coupling parameter space of T and pL q . Details of our study 
will be described in 
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